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ABSTRACT 

We determine the evolution of the faint, high-redshift, optical luminosity function of Active Galactic 
Nuclei (AGN) implied by several observationally-motivated models of the ionizing background from 3 < 
z < 5. Our results depend crucially on whether we use the total ionizing rate measured by the proximity 
effect technique or the lower determination favored by the flux decrement distribution of Lya forest lines. 
Assuming a faint-end luminosity function slope of 1.58 and the SDSS estimates of the bright-end slope 
and normalization, we find that the luminosity function must break at = —24.2, —22.3, —20.8 at 
z = 3, 4, 5 if we adopt the lower ionization rate and assume no stellar contribution to the background. The 
breaks must occur at Mg = —20.6, —18.7, —18.7 for the proximity effect estimate. Since stars may also 
contribute to the background, these values are lower limits on the break luminosity, and they brighten 
by as much as ~ 2 mag if the escape fraction of ionizing photons from high-z galaxies is consistent with 
recent estimates: / osc = 0.16. By comparing our expectations to faint AGN searches in the HDF and 
high-z galaxy fields, we find that typically-quoted proximity effect estimates of the background imply 
an over- abundance of AGN compared to the faint counts (even with / esc = 1). Even adopting the lower 
bound on proximity effect measurements, the stellar escape fraction must be high: / csc ^0.2. Conversely, 
the lower flux-decrement-derived background requires a smaller number of ionizing sources, and faint 
AGN counts are consistent with this estimate only if there is a limited stellar contribution, / csc s0.05. 
Our derived luminosity functions together with the locally-estimated black hole density suggest that the 
efficiency of converting mass to light in optically-unobscured AGN is somewhat lower than expected, 
esO.05 (all models). Comparison with similar estimates based on X-ray counts suggests that more than 
half of all AGN are obscured in the UV/optical. We also derive lower limits on typical AGN lifetimes 
and obtain ?:10 7 yr for favored cases. 



1. INTRODUCTION 

Among the long-standing goals in extragalactic astron- 
omy is to explain and characterize the population of Active 
Galactic Nuclei (AGN). Their large luminosities, compact 
sizes, and association with radio jets have lead to the as- 
sumption that AGN are powered by accretion onto super- 
massive black holes (Salpeter 1964; Zel'dovich & Novikov 
1964; Lynden-Bell 1969). Although this framework pro- 
vides a theoretical starting point, there are many questions 
that remain largely unresolved. These include explaining 
the origin of the central black holes (e.g. Eisenstein & Loeb 
1995; Madau & Rees 2001; Koushiappas, Bullock, & Dekel 
2002), understanding the fueling process, lifetime, and effi- 
ciency of the central engine (see, e.g. Rees 1984; Koratkar 
& Blaes 1999), and, ultimately, determining how quasar 
activity fits within our cosmological theory for structure 
formation. 

For many years, the role of AGN in structure formation 
was believed to be that of a tracer population, important 
in their own right, but cosmologically interesting mainly 
for their contribution to the UV ionizing background (and 
in their ability to track the collapse of structure). Recent 
indications have changed this view dramatically. It now 
seems likely that AGN play an important role in the for- 
mation of galaxies. In a reversal of sorts, this paper focuses 
on using the observed ionizing emissivity at high-redshift 
in order to constrain the evolution of the AGN luminosity 
function. Derived in this way, our luminosity functions 
relate directly to the long-standing desire to pinpoint the 
dominant ionizing sources in the Universe, and addition- 
ally help constrain models that attempt to explain AGN 
within a cosmological context. 



The AGN luminosity function has long served as a bench- 
mark for understanding the formation and evolution of 
quasi-stellar objects (QSOs) 1 (Efstathiou & Rees 1988; 
Carlberg 1990; Haehnelt & Rees 1993; Cavaliere, Perri & 
Vittorini 1997; Haiman & Loeb 1998; Richstone et al. 1998, 
Haehnelt, Natarajan, & Rees 1998; Cattaneo, Haehnelt, 
& Rees 1999; Haiman, Madau, & Loeb 1999, Kauffman & 
Haehnelt 2000; Haiman & Hui 2001; Haehnelt & Kauffman 
2001; Steed, Weinberg, & Miralda-Escude 2002). Recent 
indications that AGN activity is linked in a fundamen- 
tal way with the formation of galaxies make these studies 
all the more relevant (Heckman et al. 1984; Sanders et 
al. 1988; Kormendy & Richstone 1995; Sanders & Mirabel 
1996; Boyle & Terlevich 1998; Dickinson et al. 1998; Magor- 
rian et al. 1998; Richstone et al. 1998; Laor 1998; Wandel 
1999; van der Marel 1999; Franceschini et al. 1999; Mathur 
2000; Canalizo & Stockton 2001; Levenson, Weaver, & 
Heckman 2001; Ferrarese 2002). Specifically, the rela- 
tion between black hole mass and bulge velocity disper- 
sion (Gebhardt et al. 2000a, 2000b; Ferrarese & Merritt 
2000; Ferrarese et al. 2001) is so tight that it seems im- 
possible to understand without some significant cross-talk 
(in the form of feedback) between the AGN phase and the 
formation of the galaxy and its stellar bulge. 

Modelers attempting to understand these relations have 
been forced to test and refine their assumptions by com- 
paring to the low-redshift AGN luminosity function, or to 
bright quasar counts at high-redshift, because the faint 
population of AGN is relatively unconstrained at high-z. 
This lack of knowledge about low-luminosity objects al- 

We use the terms AGN and QSO interchangeably. In common 
parlance, a QSO is a high luminosity AGN (Mb "S —23). 
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lows considerable, unwanted freedom for model builders. 
For example, Haiman & Loeb (1998) have explored the 
idea that faint AGN are linked in a simple way with low- 
mass cold dark matter (CDM) halos. This predicts a large 
number of low-luminosity systems at high-z because less 
massive halos are relatively abundant at early times. An- 
other possibility is that AGN activity in small halos is 
suppressed by feedback processes. 2 A recent example of 
this idea is explored in Kauffman & Haehnelt (2000), who 
utilize a qualitatively plausible feedback scheme to model 
black hole properties and to match the observed evolution 
in AGN number density from z ~ to z ~ 2. Because 
their model relies on feedback that scales with host popu- 
lation and redshift, a high-z constraint on the number of 
dim objects would serve as a useful test. 

The AGN luminosity function (LF) is typically writ- 
ten as <fr(L, z)dL, and is defined as the number of objects 
per unit comoving volume at redshift z, with luminosity 
between L and L + dL. In the optical, the majority of 
studies use B magnitudes. So unless otherwise stated, L 
will denote B band luminosity throughout this paper. For 
low-redshift AGN, <f> is well represented by a broken power 
law 
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which has a break at luminosity L*, a characteristic num- 
ber density and asymptotic faint and bright slopes of 
7/ and 7b respectively. As will be discussed in §2, out to 
z ~ 2.5 the AGN LF seems to evolve only in luminosity, 
in the sense that L* gets brighter with increasing z, while 
the other parameters stay fixed (7/ ~ 1.6, 7b ~ 3.4, and 
cf)* ~ 10 3 Gpc~ 3 ). This kind of evolution is known as pure 
luminosity evolution (PLE), and the current best-estimate 
for L Sf (z) under this assumption has it rising dramati- 
cally from its local value of ~ 1O 11 L at z = to nearly 
~ 10 13 L o at z = 2.5 (Boyle, Shanks, & Peterson 1988; 
Koo & Kron 1988; Hewett et al. 1993; Pei 1995a; Boyle et 
al. 2000). Beyond this redshift, only the brightest quasars 
have been observed in significant numbers, and there is as 
of yet no evidence for a break in the luminosity function. 
In terms of the double power-law (1), the high-z obser- 
vations only measure the bright-end slope 75 ~ 2.6 — 2.9 
(Schmidt, Schneider & Gunn 1995, hereafter SSG; Fan et 
al. 2001a) and fix an overall integrated normalization that 
roughly imposes a constraint on the quantity <j)< t Lj b ~ as a 
function of z (see §2). Most interestingly, the space density 
of the observable bright quasars falls steadily from z ~ 3 
out to z - 6 (Warren, Hewett, & Osmer 1994 [WH094]; 
Kennefick, Djorgovski, & de Carvalho 1995 [KDC95]; SSG; 

2 For example, the binding energy of a galaxy of mass /jM in 
a halo of mass M and circular velocity V oc Af 1 / 3 will scale as 
£ ga l ~ 0.5f b MV 2 oc A/ 5 / 3 . Compare this to the energy released by a 
black hole of mass M, oc M shining for a time t agn at a fixed fraction 
A of its Eddington luminosity L^dd °c M,\ -E agn = AZ/£i agn oc M. 
This 

/Egal Atf 2 / 3 , suggesting that AGN feedback 
should be more important for low-mass halos. If we insert appro- 
priate numbers, we find that the energy released by a bright AGN 
over its lifetime should be comparable to the binding energy of a 
galaxy-sized halo: £ agn ~ 1.5 X 10 16 Mq km 2 s" 2 (M./10 8 Mq) 
(< agn /10 7 yr) (A/0.1), while £ gal ~ 2 X 10 15 Mq km 2 s" 2 (M/10 12 
Mq) (V/200 kms -1 ) 2 (f/0.1). Although just how this energy might 
manifest itself as a suppression mechanism is unclear, the energetics 
suggest that a significant amount of feedback is plausible. 




Fig. 1. — Schematic representation of the limits on <f>* and L* for 
2 = 4. The integrated LF (<£(< Afg)) basically constrains <f>„L2 b ~ 1 ■ 
We assume that all reasonable choices (shaded region) for the bright- 
end slope are bracketed by the values from SSG (2.87) and SDSS 
(2.58). Since the observations at high redshift show no sign of a 
break for Afg < —26, we do not consider this area of parameter 
space (vertical line). The filled square is the point in parameter 
space consistent with both the integrated LF and the parameters 
from 2dF (i.e. ~/ b = 3.41 and </>* = 1070Gpc -3 ). We also bring 
in upper limits from ionizing background measurements: proximity 
effect ("PRX.") and flux decrement ("F.D."). 



Fan et al. 2001a), but this decline cannot be faithfully rep- 
resented in terms of PLE, as the data shows 7b is flatter at 
early times. This leaves us at high redshift without a nat- 
ural extension of the bright-end LF to fainter magnitudes. 
Presumably, future QSO surveys will detect a break in the 
LF at z a3 and measure a faint-end slope. Until then, it is 
useful to examine other constraints. 

A popular technique for constraining a population of 
unresolved sources is to set an upper limit based on their 
contribution to the diffuse background light. For instance, 
AGN are strong X-ray emitters, so measurements of the 
cosmic X-ray background might be used to provide upper 
bounds on the density of AGN. The problem with this is 
that the X-ray background is measured locally, and we are 
interested in constraining a high-z population that con- 
tributes only a small fraction to the z = signal (see 
e.g. Hasinger 2002). What is preferable is a measurement 
of some background at high z by an indirect method. It 
turns out that the UV ionizing background is ideal for 
this purpose. As discussed in §3, the hydrogen ionizing 
emissivity can be measured at high-redshift by studying 
Lyman alpha absorbers along the line of sight to distant 
quasars. This measurement is especially useful because 
the derived background at a specific redshift z is roughly 
local: there is very little contribution from higher redshift 
sources because the mean free path to photo-electric ab- 
sorption is short compared to cosmological distances (see 
Madau, Haardt, & Rees 1999 [MHR] and our Appendix). 

In what follows we use this idea to construct AGN lu- 
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minosity functions that reproduce the ionizing background 
measurements, allowing in some cases for significant con- 
tributions from other sources (e.g. stars). We can demon- 
strate our general program using some simple approxima- 
tions. Let us assume (as we do throughout) that we can 
approximate each AGN as emitting light with a self-similar 
spectrum that varies only in normalization from object to 
object (§2.2). This implies a fixed ratio between an AGN's 
specific luminosity at the 912A Lyman edge to its luminos- 
ity in the B-band: n u = L v /L ~ 10 18 erg/s/Hz/L Q . The 
comoving emissivity at the Lyman edge is then: e„ = n v l, 
where I is the luminosity density of AGN (in the B-band). 
For a reasonable LF, I will be dominated by objects near 
the break, so we expect £ ~ 4>*L*, or e v ~ rj u <j)^L^ . A 
slightly more accurate estimate of I comes from integrat- 
ing the LF (1) from L = to oo, and with this we obtain 
e v ~ r/^0*i*[(2 - 7/) _1 + (76 - 2) -1 ]. Thus, a measure- 
ment of the ionizing emissivity at some redshift mostly 
constrains the parameter combination </>*L*, with a weak 
dependence on the LF slopes (as long as neither slope ap- 
proaches 2). Direct observations, on the other hand, mea- 
sure the bright-end slope jj, and fix the normalization of 
the bright-end tail. Together, then, the limits from e v and 
direct observations can effectively restrict the acceptable 
values of and L*. 

A schematic illustration of this is shown in Figure 1, 
where we plot the available parameter space at z = 4. Di- 
rect observations of the density of bright QSOs leads to 
a constraint of cj>„L* oc Lt~~' b (see §2.1), which we draw 
(diagonal lines) for two values of 7^ from separate surveys 
(2.87 from SSG; 2.58 from Fan et al. 2001a). We presume 
that the best estimate of the bright-end slope lies between 
these two measurements, and therefore, that the values of 
<j>* and L* are situated somewhere in the shaded region. 
The fact that the break has not yet been observed down to 
the limiting magnitude Ly lm ~ 10 Lq of SDSS allows us 
to set an upper limit on L* (vertical line). We can further 
narrow down the parameter space with two different (con- 
flicting) measurements of the ionizing emissivity at the 
same redshift (the higher line comes from the proximity 
effect and the lower line comes from the flux decrement 
distribution, as discussed in §3). We draw these as upper 
limits (horizontal bands), since contributions to the back- 
ground from non-AGN sources will require fewer AGN, 
further restricting the allowed range of <^*L*. The width 
of the bands reflects the slight dependence of the emissiv- 
ity on 7b. As for the faint-end slope, we have fixed it to 
the low redshift value of 1.58, but we explore other values 
of 7/ in §5.1. What is evident in the figure is the fact that 
the emissivity measurements provide an extremely useful 
limit on the parameter space left available from the direct 
observations of bright QSOs. 

In the next section we review what is known about the 
optical luminosity function from direct observations, and 
discuss some possible caveats associated with dust and 
gravitational lensing. We also discuss our assumed tem- 
plate AGN spectrum. In §3 we review the different ioniz- 
ing background measurements, and in §4, we describe our 
calculation of the AGN contribution to the background as 
well as the contribution from stars and IGM reemission. 
We present our results in §5 for three separate models of 
AGN emissivity, and we compare the derived LFs to past 



and future faint surveys. In §6 we discuss our results in 
the context of other constraints on the relative AGN and 
stellar emissivities at high-z. We also examine what our 
results imply for the AGN efficiencies and lifetimes. Our 
conclusions are summarized in §7 

Unless otherwise stated, we assume that the cosmology 
is one with the universe made flat by a cosmological con- 
stant Q m = 1 — Cl A = 0.4, and that the Hubble parameter 
at z = is H = 65 km s' 1 Mpc -1 . 

2. AGN PROPERTIES: COMPOSITE SPECTRA AND 
OBSERVED LUMINOSITY FUNCTION 

In this section we will discuss the empirical AGN lumi- 
nosity function (LF) and go on to present our assumed 
template spectral energy distribution, which allows us to 
calculate the ionizing emissivity for any given <p(L, z). As 
mentioned previously, L will always have units B-band so- 
lar luminosities. We assume that the spectrum of an AGN 
varies only in normalization from object to object, and 
we present our assumed distribution in terms of specific 
luminosity L v in units of erg/s/Hz. 

2.1. Empirical AGN Luminosity Function 

In §1 we discussed how low- redshift observations indi- 
cate that the AGN luminosity function is well-described 
by a double power law shape given by Eq. 1. Although 
yet to be confirmed by observations, we will assume that 
this shape provides a useful characterization of the LF at 
all redshifts. Under this premise, the AGN LF has four 
main parameters that must be described at each z: 7/, 76, 
<f> t , and L*. 

Detecting AGN at low redshift is done primarily by se- 
lecting for strong UV emitters (e.g. Hartwick & Schadc 
1990; Hewett et al. 1995). The largest survey to date is 
that of the 2dF (Boyle et al. 2000), which tracks the evo- 
lution of QSOs over the redshift range 0.35 < z < 2.3. 
Their data is consistent with the assumption of pure lumi- 
nosity evolution (PLE). Under PLE, all of the LF evolu- 
tion is contained in L*, while the shape (7/, 7;,) and nor- 
malization (</>*) remain constant. Their fit for a ACDM 
cosmology has 7/ = 1.58, 76 = 3.41, 0* = 1070Gpc~ 3 , 
and log 10 (i*) = 11.24 + 1.36z - 0.27z 2 . The rapid rise 
in L*(z) characterizes an increase in the number density 
of AGN as we look to higher z. We plot the 2dF fit in 
Figure 2 for several redshifts. For reference, we have only 
plotted the fits over the range of (absolute) luminosities 
probed by the 2dF. The rapid evolution and characteristic 
LF shape found by the 2dF for zs2.5 is consistent with 
previous work in overlapping redshift ranges (see e.g. Pei 
1995a). 

Beyond z <~ 2.5, however, our knowledge about the AGN 
population is much less complete. The information we do 
have is based on the bright end of the LF, which is ob- 
served to fall off gradually out to z>5 (WH094, SSG, Fan 
et al. 2001a, b). Of course, if PLE holds to high-z then 
the evolution of the brightest quasars is sufficient to de- 
fine the entire LF evolution. For example, under the PLE 
assumption, Pei (1995a) and later MHR used the observed 
bright-end evolution (along with low-z data compilations) 
to extrapolate the evolution of the AGN LF out to high- 
z. Models of this kind have been popular for estimating 
the number of high-redshift AGN and calculating expected 
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Fig. 2. — Fits to the observed LF for several redshifts. The 
three lowest redshift measurements come from Boyle et al. (2000) for 
16.50 < M B < 20.85. The two highest come from Fan et al. (2001a) 
for 18 < i* < 20. 



AGN contribution to the ionizing background. Yet even 
before the most recent SDSS data, there were indications 
that the PLE assumption might break down at early times. 
Specifically, the bright-end slope obtained by SSG (and by 
KDC95 at slightly larger magnitudes) is flatter (7;, ~ 2.9) 
than that measured for local AGN. The recent SDSS data 
(Fan et al. 200fa) reveal an even flatter bright-end slope 
(7b ~ 2.6) for 3.5 %z §5.0. And although not as statistically 
significant, the data at slightly fainter magnitudes from 
the Isaac Newton Telescope Wide Angle Survey (Sharp et 
al. 2001) are consistent with the results from SDSS only 
if 7b s2.9. All together, these observations provide strong 
evidence that the LF evolves in shape as well as in lumi- 
nosity, and therefore PLE does not extend to high redshift. 
In this case, constraints on the faint-end of the LF become 
all the more valuable. 

We will normalize all of our high- z LF models to match 
the SDSS results (Fan et al. 2001a, b) over the range in 
L that they are measured. The best-fit LF of Fan et 
al. is shown in Figure 2 for z — 3.6 and 5.0 over the rele- 
vant range of absolute magnitudes. The best-fit bright-end 
slope is 7b = 2.58 (the quoted error is ±0.23) and the num- 
ber density evolution over the redshift range 3.5 £z §6.0 
follows the integrated constraint 



$(M B < -26, z) = lo 1 "" - 47 ^" 3 ) Gpc 



-3 



where 



$(< M B ,z) 



L{M E 



4>(L, z)dL, 



(2) 



(3) 



and log L (Mb) = -0A(M B -M 0)B ). Note that for 7 b > 1 
and L(Mb) 2> the integral (3) approximates to ~ 
</)*(L*/L m ) 7i>_1 /(7ft — 1). The parameterized normaliza- 
tion from the SDSS work is shown in Fig. 3 along with data 
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Fig. 3. — The number density of bright quasars for z > 3. The 
solid line is one of our models constrained to fit the SDSS data, and 
the dot-dash line is the PLE LF from MHR. The data points are 
from Hartwick & Schade (1990), Warren, Hewett, Osmer (1994), 
Schmidt, Schneider, Gunn (1995) and SDSS (Fan et al. 2001a). 
Since all the of the compiled data sets assume Qj^ = 1, we have 
scaled our fit to match this cosmology. Note also that for Hartwick 
& Schade and Warren, Hewett, Osmer the error bars are summed 
over the differential LF uncertainties and, therefore, are likely over- 
estimated. 



from several other high redshift surveys. Except for the 
seemingly anomalous point from WH094 at z = 3.2, all 
the integrated LF measurements appear to agree. The fig- 
ure also illustrates how the PLE model of MHR disagrees 
with the recent data, and highlights the need to recon- 
sider the AGN contribution to the ionizing background at 
high-z. In what follows we will use the SDSS result to fix 
the bright-end slope and normalization of our constrained 
AGN LFs. 

Before going on, we should address possible selection 
effects which could affect the empirical LF at high red- 
shift. One possibility is that dust, both intrinsic to the 
host galaxy and in the foreground, may be affecting the 
observed evolution as well as shape of the AGN LF. For 
example, Fall & Pei (1993) showed that damped Lya sys- 
tems (DLAs) could be blocking 10-70% of the bright op- 
tical QSOs at z ~ 3. But the CORALS survey (Ellison 
et al. 2001) found that the distribution of DLAs was ba- 
sically the same in both radio and optically selected QSO 
samples. Since dust should not affect radio wavelengths, 
this implies that intervening dust is not a significant bias. 
Moreover, Fan et al. (2001a) find the distribution of spec- 
tral indices in the SDSS sample is similar to that at low 
redshift, which would be unlikely if there was a substantial 
amount of dust along the line of sight. So too, any increase 
in reddening with redshift would presumably be at odds 
with the close agreement in LF evolution from surveys 
with very different selection techniques: the broad-band 
color selection in SDSS; the Lya-emission selection of SSG; 
and even radio selection (e.g. Hook, Shaver & McMahon 
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1998), which should be representative of full sample, as- 
suming that the ratio of radio-loud to radio-quiet QSOs is 
redshift independent (Stern et al. 2000). However, none 
of this would seem to preclude the existence of a distinct 
population of highly extinguished QSOs, only observable 
in the X-ray (and perhaps the IR). We come back to this 
possibility when we address the relic black hole density in 
§6.2. 

Another effect that could distort the LF is gravitational 
lensing (see Blandford & Narayan 1992 and references therein). 
Because of the shape of the LF, the magnification bias is 
stronger for more luminous objects, i.e. it is more likely 
that the bright-end of the LF is contaminated by arti- 
ficially brightened AGN. According to Pei (1995b), the 
number of objects observed with L = 10L* that are lensed 
could be as much as 44% at z — 3 and 68% at z = 5. 
However, these large lens fractions come from models that 
assume a now unreasonable density of compact objects 
(on the order of 10% the critical density). More recent 
estimates of the fraction of lensed objects in magnitude 
limited surveys are less than ~ 20% (Barkana & Loeb 
2000; Wyithe & Loeb 2002). Because these are relatively 
small corrections, we will ignore the effect of magnifica- 
tion bias in our bright-end normalizations. However, it 
should be pointed out that the lens fraction gets larger the 
farther an observed luminosity is above the break. That 
is, the probability that the QSOs observed in the SDSS 
(L > 10 12 5 L Qi b) are lensed increases the smaller L* turns 
out to be. In our emissivity constrained models (see Figure 
1), we explore break luminosities that are several orders of 
magnitude below those investigated by Pei (1995b). How- 
ever, these values of L* are probably unrealistic because 
they seem to be in conflict with some of the faint AGN 
searches (see §5). 

2.2. Lyman Limit emissivity and Composite Spectra 

Given an optical AGN LF 4>(L, z), we are concerned with 
calculating the implied emissivity of ionizing photons. If 
e(v, z) is the comoving emissivity of photons at frequency 
v and redshift z, then we can write 

/>oo 

e(v,z) = J dL-<j>(L,z)-L v (L,v). (4) 

Here (and throughout this work) we adopt a minimum 
AGN luminosity L min = 1O 8 L , B (or M B = -14.5). This 
choice is comparable to the nuclear luminosity of the faintest 
Seyferts (Londish et al. 2000). Fortunately, the choice 
of L m i n does not strongly affect our results as long as 
L m in L*. 3 Expression 4 provides a link between the 
emissivity and <j){L,z), once we specify the input spec- 
trum, £„(L, v). We will assume that, on average, the spe- 
cific luminosity of an AGN (as a function of frequency v) 
varies only in normalization such that L v oc L, with no 
z— dependence (see Kuhn et al. 2001). 

From the UV to optical, an AGN spectrum, L v , is rea- 
sonably well approximated by a double power law with 
a break near Lya at 1216A. Since the Galaxy is opaque 
to UV photons, we must rely on high redshift quasars to 

3 If we let L m ; n = 10 9 Lq b we find it affects our results (§5) by less 
than 30%. And making the minimum smaller than our fiducial value 
has almost no effect, as the integrals converge. 
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Fig. 4. — Our assumed template AGN spectrum shown along 
with the composite spectrum of Elvis et al. (f994). For reference, 
we have drawn vertical lines at the Lyman limit and a shaded band 
to indicate the wavelength range of the B-band. 



determine a typical quasar spectrum short-ward of the Ly- 
man limit. For A < 1200A, we will use results from the 
most recent UV survey from HST (Telfer et al. 2002). 
Telfer and collaborators find that the specific luminos- 
ity of (radio quiet) AGN scales as L v oc v~ auv , ajjy — 
1.57 ± 0.17, from A - 500 to - 1200A. 4 Long-ward of 
Lya, we use the result of Vanden Berk et al. (2001), who 
relied on a composite of quasar spectra from the SDSS to 
determine L v oc z^ 44 from 1200A to - 5000A. 

Figure 4 shows our assumed spectrum compared to the 
composite spectrum of Elvis et al. (1994). The specific 
luminosity is normalized relative to the B-band luminosity 
L. For our B-band transformation we follow SSG and use 
M B = M A b (4400A) +0.12. Of course, what we are mainly 
interested in is the ratio of L v at 912A [y = vu) to the 
optical luminosity L, which for our spectrum is 

s M^h) = 10 i8. 05 ergs/s/Hz 
L L 0! b 

From combining the optical and UV power laws, we esti- 
mate an uncertainty in r\ v of ±0.10 in the exponent. This 
level of variation will not strongly affect our conclusions. 
For reference, the value implied by Elvis et al. (1994) is 
rj v — 10 1815 in the same units, and the value obtained 
by Shull et al. (1999) using a sample of 27 Seyferts is 
r\ v = 10 17 93±0 04 . Note that Shull et al. also find a slight 
luminosity dependence on this ratio, but we will ignore 
that possibility here. With our assumed L v spectrum, Ex- 
pression 4 fully describes the relationship between a given 

4 We assume that the minority population of radio loud objects will 
contribute very little to background. It is also worth noting that 
the Telfer et al. (2002) slope is slightly harder than previous results 
from HST (Zheng et al. 1997), which found a vv = 1.83 ± 0.15. If 
we rerun our analysis with this steeper continuum, we find that our 
results change by less than 20%. 
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optical LF and the ionizing emissivity as a function of red- 
shift. 

3. IONIZING BACKGROUND MEASUREMENTS 

From the Gunn-Peterson effect (Gunn & Peterson 1965), 
it is clear that the IGM is highly ionized. If it can be 
assumed that this high ionization state is due to photo- 
ionization (as opposed to collisional excitation), then mea- 
surements of the IGM, particularly the Lya forest, can re- 
veal information on the nature of the ionizing background. 
We can characterize this background with the hydrogen 
ionization rate: 



\z) = 4tt f 

J V 



50 , J{y,z) 
dv an(y). 

1>YL 



V 



(6) 



Here the cross section to hydrogen photo-electric absorp- 
tion is <th = 6.35 x 10~ 18 cm 2 • (^/^h)~ 3 and J(y,z) is 
the background intensity (in units of ergs/s/cm 2 /Hz/sr). 
As will be discussed in the next section, the ionizing in- 
tensity relates simply to the AGN (and stellar) emissivity, 
J(z) oc s(z), so we may approximate the spectrum of back- 
ground intensity as J oc v~ aiJV . With this approximation 
we obtain: 



r_i 2 (s) = ™ J-2l(z) 

6 + Oiuv 



(7) 



where J-2i(z) = J(vk, z) in units of 10 _21 ergs/s/cm 2 /Hz/sr 
and r_i2 = r/10 _12 s _1 . If there were only a single domi- 
nant source to the background (e.g. QSOs), then we could 
take auv to be equal to the typical spectral index of that 
source population. We will however be exploring cases 
where there are more than one dominant type of emitters, 
and so we will not immediately assume a value of auv- 5 
For this reason, in the next section, we will quote mea- 
surements of the ionizing background using r_i 2 (instead 

Of J-2l). 

Under the assumption of photo-ionization equilibrium, 
r should scale as the ratio of ionized to neutral hydrogen 
(nmi/n-m)- Given this ratio one may infer a value for T. 
At high-z, two techniques exist for measuring the amount 
of neutral hydrogen, and they both rely on using Lya lines 
in the spectra of distant QSOs. 6 One technique relies on 
analytic modeling of the expected distribution of lines in 
the vicinity of a QSO (the proximity effect) and the other 
uses hydro-dynamical simulations to model the expected 
distribution of lines along the line of sight to the quasar 
(what we will call the "flux decrement" analysis). We use 
the following two subsections to summarize the results of 
each technique. 

3.1. Proximity effect 

The proximity effect was first discussed by Bajtlik, Dun- 
can, & Ostriker (1988). The technique relies on the ob- 
served decrease in the number of Lya lines in the vicinity 

5 Moreover, the processing of the background through absorption 
and reemission in the IGM will alter its spectral shape (Miralda- 
Escude & Ostriker 1990). 

6 Note, because the Lya forest is very thin for z < 1.7, it is difficult 
to make measurements of the ionizing background in this way at low 
rcdshift (e.g. Shull et al. [1999] and references therein). Although 
see Dave & Tripp (2001). 



of a QSO relative to what one would have expected in the 
absence of the QSO. In principle, this decrease is a result 
of ionizing radiation from the QSO itself, which tends to 
reduce the neutral fraction in nearby absorbers. The num- 
ber of observable lines should decrease correspondingly. 
An estimate of the ionizing rate can be obtained by deter- 
mining the distance from the quasar at which the number 
of lines is equal to the background expectation. Over the 
years, measurements of J-21 (assuming auv ~ 1-8) have 
taken values between ~ 0.7 and ~ 3 (Williger et al. 1994; 
Bechtold 1994; Fernandez-Soto et al. 1995; Cristiani et 
al. 1995; Giallongo et al. 1996; Cooke, Espey & Carswell 
1997; Scott et al. 2000; Liske & Williger 2001). Since what 
is effectively measured is a ratio of proper distances, this 
analysis should be relatively independent of cosmological 
parameters. Most estimates of the proximity effect assume 
that the background is constant over the redshifts being 
sampled, so this technique does not appear to uniquely de- 
termine the evolution of r_i 2 (2;). Haardt & Madau (1996) 
and Fardal et al. (1998) have tried to fit the whole set of 
proximity effect measurements with a Gaussian. Scott et 
al. (2000) claim that their observations are well fit by the 
following parameterization from Fardal et al. (1998): 



r_ 12 (z) = 1.2(1 + z) a58 exp 



-(z-2.77) 2 
2.38 



(8) 



This fit is represented by the dot-dashed line in Figure 5. 

There is some concern that the value of T_ 12 determined 
from the proximity effect is overestimated. This is be- 
cause QSOs are likely to be found in environments that 
are denser than average (Pascarcllc et al. 2001; Ellison et 
al. 2002), so that the regions of excess ionization end up 
being smaller than they would have been if the region was 
of average density. Such overdensities would likely scale 
with the mass/luminosity of the AGN, which might ex- 
plain the slight anti-correlation between more luminous 
QSOs and their proximity effects (see Cooke, Espey & 
Carswell 1997). Loeb & Eisenstein (1995) claim this bias 
could cause the proximity effect measurements to overes- 
timate the background by as much as a factor of 3. The 
flux decrement technique discussed in the next subsection 
does in fact yield a lower ionizing background than the 
proximity effect, as would be expected if these concerns 
are valid. However, we will keep an open mind about the 
issue and discuss the constraints on faint AGN associated 
with each technique separately. 

3.2. Flux Decrement 

Another way to utilize the Lya forest to measure the 
ionizing background requires making a theoretical predic- 
tion as to the amount of (unobserved) ionized hydrogen in 
each "cloud" along the line of sight to a quasar. This can 
be done using CDM theory (and N-body codes) coupled 
with hydrodynamical simulations that model gas evolution 
(Cen et al. 1994; Zhang et al. 1995; Hernquist et al. 1996; 
Miralda-Escude et al. 1996; Rauch et al. 1997; Dave et 
al. 1997; Wadsley & Bond 1996; Zhang et al. 1998; Mc- 
Donald et al. 2000). These simulations essentially gener- 
ate a distribution of hydrogen Lya optical depths which 
can be mapped to a distribution of flux decrements (F = 
exp(— t)). The distributions generated in this way con- 
vincingly reproduce observations, adding to the long list 
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Fig. 5. — The ionizing rate as a function of redshift. The solid dots 
are from hydro-simulations by McDonald & Miralda-Escude(2001). 
A linear fit has been drawn through the points. The solid boxes and 
open square are representative of recent proximity effect measure- 
ments (see Scott et al. 2000), where we have used the conversion 
T_i2 = 2.64 ■ J_2i from Fardal et al. (1998). We have also drawn 
the Gaussian fit to this data from Fardal et al. 



of successes for the CDM paradigm. Although the shape 
of the flux decrement distribution depends upon details 
of the cosmology as well as models for the temperature- 
density relation, the normalization is proportional only to 
the baryon density, the ionization rate 7 , and the Hubble 
parameter: 



{n b h 2 f 
T{z) H( z y 



(9) 



Often times this relation is used to put constraints on 
fl b h 2 (e.g. Hui et al. 2002). But McDonald & Miralda- 
Escude (2000) [hereafter M&M-E] assume the baryon frac- 
tion from BBN: n b h 2 = 0.02 (Walker et al. 1991; Buries & 
Tytler 1998), so as to measure T(z) from a sample of QSO 
spectra. In a ACDM cosmology, over the redshift range of 
~ 2.3 — 5, their results can be approximated by 



-12 



_ i n-0.242+0.4 



(z) = 10 



(10) 



We plot this flux-decrement result along with the proxim- 
ity effect results in Figure 5. The flux decrement determi- 
nation is a factor of a 4 below that of the proximity effect. 
We discuss how these differences affect constraints on the 
AGN LF in §5. 8 

7 It is usually assumed that the ionization rate is homogeneous. 
Gnedin & Hamilton (2002) included a spatially inhomogeneous ra- 
diation background from galaxies in their simulation and found that 
the mean background came out larger. This could mean that the 
flux decrement measurements underestimate the ionization rate by 
at least 20%; however, a significant contribution from AGN would 
likely lessen this effect. 

8 In recent papers (Cen & McDonald 2002; Fan et al. 2002), the flux- 
decrement analysis has been extended to z ~ 6, and evidence for the 
epoch of reionization has been discussed. If the epoch of reionization 



4. THE IGM: AGN AND STARS AS IONIZING SOURCES 

Now that we have measurements of the background ion- 
ization rate, we need to relate them to the LF. This rela- 
tion comes about via Eq. 7 which connects T to the ionizing 
intensity. In principle, the intensity at any time is made 
up of the integrated contribution of all sources over the 
history of the universe. However, as discussed in detail by 
MHR, the contribution from distant sources is significantly 
degraded by attenuation. The intensity can be written as 
an integral over the (comoving) emissivity of sources as a 
function of redshift: 



J{y,z) 



(1- 



dz—-£(D,z)-e- T ° a <- u < z > s \ (11) 
dz 



where v = v (l-\-z)/(l + z), and r e g is the effective optical 
depth due to photo-electric absorption in the IGM. In the 
Appendix we discuss how we model the IGM and r e / / fol- 
lowing the prescription of MHR. In practice we calculate 
J using the full integral expression, but for purposes of 
illustration, it is useful to make the following approxima- 
tion. If we assume that only ionizing sources within one 
absorption length, Al (where r e ff(AZ) = 1), contribute to 
the background, then Eq. 11 reduces to Jiy.z) ~ (1 + 
z) 3 £v(z)Al/4ir. At the Lyman edge, our adopted IGM 
model (Appendix) gives Al ~ 39Mpc[(l + z)/4]~ 4 - 5 . Since 
Al <C c/H(z) this is often called the "local source" ap- 
proximation, and it yields: 



J_ 2 i(z) ~ 0.17 (1 + z)- 1 - 5 E2i{v H , z). 



(12) 



Here we have introduced the symbol £24 which is the emis- 
sivity in units of 10 24 ergs/s/Hz/Mpc 3 . 

Equations 11 and 7 allow us to relate an observed ion- 
ization rate T(z) to a background ionizing emissivity e{z). 
However, the relation depends sensitively on the spectral 
slope of the background: Y oc J/(3 + auv) z/{3 + auv)- 
Because the far- UV slopes of galaxies and AGN are signif- 
icantly different, deriving a limit on the background emis- 
sivity based on the ionization rate would necessitate as- 
suming something about the background population. We 
would prefer a more general parameter that allows us to 
consider stars as well as AGN as major contributors to 
the background. For this purpose we propose the follow- 
ing "weighted" emissivity parameter 



; 24 



(13) 



l uv 



which facilitates a fair comparison between different pop- 
ulations, i. This definition allows us to rewrite the (ap- 
proximate) Equation 12 as 



r_ 12 (z)^2.0 (l + z)- 1 - 5 £ 



(14) 



where the sum is over all populations of ionizing ionizing 
sources, i = AGN, stars, reemission, etc. This expression 
can be rearranged to place an upper limit on the weighted 

has been detected, then the drop in the ionization rate at z ~ 6 will 
be due in part to a decrease in the mean free path of continuum 
photons, and not merely to a change in the ionizing emissivity. See 
§4. 
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emissivity coming from any individual population. In par- 
ticular, we are concerned with limiting the emissivity from 
AGN: 

e Q (z) ~ 0.5r_ 12 (z) (1 + z) 1 - 5 - [£*(z) + e R ...]. (15) 

The bracketed terms on the right hand side are included 
to illustrate how the allowed AGN emissivity is reduced by 
the presence of additional contributing populations. For 
example, i* represents the weighted emissivity from stars 
(see §4.1) and e R represents the contribution from cloud 
reemission (see §4.2). The contribution from clusters of 
galaxies is likely negligible (Randall & Sarazin 2000). 

Fig. 6 illustrates the upper limit on weighted AGN emis- 
sivity i imposed by the proximity effect (short dashed 
line) and flux decrement (long-dashed line) measurements 
of T(z). As will be the case throughout, we have used 
the full integral expression given by Eq. 11 in order to 
place limits on the weighted emissivity (as opposed to the 
approximate expressions used for illustrative purposes in 
Eqs. 12, 14, and 15). However, in the language of Eq. 15, 
these two constraints amount to setting all of the brack- 
eted (non-AGN) emissivities to zero. The thin solid line 
(Model A), on the other hand illustrates how the flux- 
decrement-derived limit on the AGN emissivity changes if 
we include a (conservative) estimate of the stellar emissiv- 
ity, e*. We describe this estimate in the next subsection. 
For reference we show the contribution to § arising from 
AGN that have been directly observed by the 2dF at low- 
z (bold solid line), and by SSG and the SDSS at high-z 
(data points). 

We now go on to discuss our stellar emissivity estimate 
in more detail. 

4.1. Emission from Galaxies 

In addition to AGN, star-forming galaxies are an ob- 
vious contributor to the ionizing continuum. At energies 
above the Lyman edge, a star-forming galaxy spectrum 
is dominated by hot, short-lived O stars (to < 10 7 yr). 
Theoretically, the number of Lyman continuum photons 
emitted by an (unobscured) L* galaxy should be about 
10 -3 smaller than the number emitted from an AGN 
(Madau & Shull 1996). However, the space density of L* 
galaxies is a factor of ~ 10 4 higher than that of L* AGN 
at z ~ 0, so stars could conceivably dominate the ionizing 
background. What mitigates this calculation is the frac- 
tion of Lyman continuum photons that can actually escape 
from the galaxy, / csc . 

Theoretically, cold gas in a galaxy could trap most of 
the Lyman continuum emission (e.g. Haehnelt et al. 2001 
and references therein). In agreement with this expecta- 
tion, most local searches for Lyman continuum emission 
from galaxies have come up empty. At low redshift, using 
the Hopkins Ultraviolet Telescope, Leitherer et al. (1995) 
found that nearby starburst galaxies had a very small es- 
cape fraction, / osc < 3% (see also Hurwitz et al. 1997). 
Similar upper limits can be found from FUSE (Deharveng 
et al. 2001) and from HST measurements in the HDF (Fer- 
guson 2001). The escape fraction from the Milky Way can 
be estimated from Ha measurements in the Magellanic 
Stream (/ csc ~ 6%), but there is some uncertainty associ- 
ated with this estimate (see Bland-Hawthorn & Maloney 
2001). 




1 2 3 4 5 

redshift 

Fig. 6. — The weighted emissivity: e(^jj)/(3 + cijyy-). The 
bold solid line is the weighted emissivity implied by observed QSOs 
(-26 < M B < -23) from 2dF (Boyle et al. 2000). The crosses are 
from SSG and the triangles are from the SDSS (Fan et al. 2001a), 
both involving a conversion from Q(Mb < —26). The star is the 
weighted LBG emissivity from Steidel et al. (2001). Limits have 
been drawn from the two ionization rate measurements: proximity 
effect (PRX.) and flux decrement (F.D.), see §3.1 and §3.2. We 
have also included our Model A (solid line), which is the allowed 
AGN emissivity accounting for a stellar contribution to the flux- 
decrement-derived background (see §5). 

At high redshift there are indications that the escape 
fraction is much higher than these local observations sug- 
gest. Steidel et al. (2001) reported the first evidence for 
galactic emission beyond the Lyman edge using a set of 
galaxies at z ~ 3. Their result is based on a sample 
~ 1000 spectroscopically-confirmcd, star-forming, Lyman 
Break Galaxies (LBGs), for which the total emissivity at 
1500A is fairly-well determined: £^ 4 (1500A) = 180 at 
z = 3 for our ACDM cosmology (Steidel et al. 1999). Us- 
ing 29 galaxies effectively taken from the bluest quartile of 
their LBGs, Steidel et al. (2001) derived a composite spec- 
trum, and used it to study the flux density at 1500A com- 
pared to 900A. They obtained the surprisingly low value 
f [1500] /f [900] = 4.6 ± 0.1, which is nearly consistent with 
theoretical expectations in the absence of any internal ab- 
sorption (~ 3 — 5 according to Leitherer et al. 1999, Steidel 
et al. 2001, and references therein). If one assumes the in- 
trinsic ratio is 3, then the escape fraction implied by their 
measurement is / GS c = 65%. If all LBGs were emitting 
ionizing photons at this rate then they would likely dom- 
inate the AGN-contribution to the background, and, in 
fact, over-produce the ionizing rate implied by flux decre- 
ment measurements. However, Giallongo et al. (2002) re- 
cently looked at 2 random LBGs and found no Lyman 
continuum photons. It seems more likely, or at least more 
conservative, to assume that only the bluest quartile of the 
LBG population at z ~ 3 is emitting ionizing photons with 
/csc =65%, and the rest are much more opaque. Under 
the assumption that only the bluest fourth have escaping 
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Lyman continuum photons, the average escape fraction for 
the entire population would be / esc = 65%/4 ~ 16%. This 
is what we will assume here. 

In terms of the escape fraction, the implied z = 3 emis- 
sivity at the Lyman edge is s^^vh) = (/esc/3) -£24 (l^OOA) c: 
60 • / csc . However, what we are concerned with is the stel- 
lar contribution to the ionizing rate, which depends on the 
far-UV background slope auv- As in Eq. 13, we charac- 
terize the slope-dependence using the weighted emissivity 
of stars: 



£24^) _ 60 -/c; 



a 



(16) 



uv 



l uv 



where we have implicitly assumed z — 3. For the slope 
we will again rely on the results of Steidel et al. (2001). 
Between 1100 A and 900A, their spectrum yields a slope 
of ajjv — 6-7, and we will assume that this power-law can 
extrapolated into the Lyman continuum without signifi- 
cant error. 9 With auv — 6.7 and / csc = 16%, the implied 
weighted emissivity for stars at z — 3 is e* ~ 1.0. This 
value is represented by the five-pointed star in Figure 6. 

One concern is that the value of auv we derived us- 
ing the Steidel et al. (2001) composite spectrum is much 
softer than theoretical models often assume. A more com- 
mon assumption is auv — 2 (e.g. Miralda-Escude & 
Ostriker 1990; Madau & Shull 1996; Shull et al. 1999; 
Hachnelt et al. 2001). While our use of the Steidel slope 
is well-motivated, if one were inclined to believe a differ- 
ent value, we could simply absorb the change in auv by 
re-interpreting our assumed / eS c- For example, if we adopt 
the theoretically-motivated value auv — 2.0, then we ob- 
tain the same e* using / esc ~ 8%. Since our main goal 
in using this stellar emissivity model is to illustrate how 
it will affect the implied AGN LF, any reader unhappy 
with our adopted slope can simply regard our assumption 
to be £24 — 1 at z = 3, and interpret / OS c accordingly. 
Nonetheless, we feel that the Steidel et al. result provides 
an observationally-motivated choice and we adopt the im- 
plied normalization for our model constraints in §5. 

In order to extend our AGN analysis to zz3, we must 
extrapolate the stellar contribution to higher rcdshift. We 
do so by assuming that the stellar emissivity evolves in 
the same manner as the star formation rate density of the 
universe: e*(z) = e*(z = 3) • p*(z) / p\(z = 3). For za2, we 
assume that the star formation rate evolves as 



p*{z) oc 



4z 



1 + (z/4.1) 4 -i 



1 



(17) 



The shape of this function provides a good fit to the star 
formation predictions of Somervillc, Primack, & Faber 
(2001) (R. Somerville, private communication), and also 

9 One might expect a particular starburst spectrum to show a break 
at the Lyman limit, so that our choice of extending the Far UV 
(A > 912A) spectral index into the Lyman continuum may be un- 
warranted. However, since, within our picture at least, the main 
contributors have almost no attenuation such an effect would not 
be expected. For example, it is possible that the bluest LBGs are 
being seen through "supcrbubblcs" (Dawson et al. 2002) in which 
the neutral hydrogen has been partially blown out or ionized away. 
This would leave the spectrum relatively smooth across the Lyman 
break, whereas along some other line of sight, the same starburst 
may have almost no Lyman continuum emission. In such a scenario, 
/esc would absorb these viewing-anglc effects. 



provides a good representation of the data, although the 
scatter in the data is large (see, e.g. Steidel et al. 1999, Poli 
et al. 2001, and the compilation in figure 12 of Springel & 
Hernquist 2002). 

4.2. Reemission from Clouds 

Before going on to the results, we mention that an- 
other likely contributor to the ionizing flux arises from the 
reprocessing by the IGM of the primary radiation from 
stars and AGN (see Haardt and Madau 1996; Fardal et 
al. 1998). This effectively transfers some high energy pho- 
tons to lower energies where they are more likely to ionize, 
in our case, hydrogen. The H 1 recombinations are right at 
the Lyman edge, so these photons are quickly redshiftcd 
to below threshold. However, He 11 recombinations and 
two photon de-excitations will have a larger contribution 
to the ionizing background (Shull et al. 1999). Since stars 
are not thought to emit many helium ionizing photons, 
this reemission will mostly come from reprocessing AGN 
radiation. Haardt & Madau found that this reemission 
increases the hydrogen ionization rate due to QSOs by 
about 40%, independent of redshift. In this case we can 
write e R = 0Ae Q . According to Fardal et al. (1998), this 
fraction is ~ 25%. For our fiducial models explored in 
the next section, we will assume that reemission is negli- 
gible (e R — 0), but in §5.1 we explore how including 40% 
reemission will affect our derived AGN LFs. 

5. RESULTS 

We will present our constraints under assumption that 
the LF takes the form of Eq. 1, and thus has four free pa- 
rameters to be constrained at each redshift: 7/, 7^, and 
L„. By forcing the bright-end to match the results of di- 
rect observations we have two constraints: the SDSS slope 
(7b = 2.58) and normalization (Eqs. 2 and 3). A third 
constraint comes from the ionizing background, which ef- 
fectively limits the integrated faint-end emissivity. We fix 
the final parameter by first setting the faint-end slope at 
its low-redshift value, 7/ = 1.58, and then we explore how 
our results change for other values of 7/ in §5.1. With 7^ 
and 7/ fixed, our emissivity constraints can be expressed 
by simply quoting the implied break luminosity L* (or the 
break magnitude Mg), since the SDSS normalization fixes 
(/>» for a given L„. 

Table 1 summarizes the ionizing background models we 
have used to derive our LF's. In Model A, we assume 
that the total ionizing rate at high-z is set by that quoted 
by M&M-E (using the flux decrement distribution tech- 
nique). Model A also assumes that starlight contributes 
to the background at a level consistent with what Stei- 
del et al. (2001) found, / osc = 0.16. Our second example 
(Model B) also assumes the M&M-E ionizing rate but now 
with a negligible stellar contribution / osc = 0.0. Finally, 
in Model C, we assume the (higher) ionizing rate inferred 
from the proximity effect analyses with / csc = 0.0 (we ex- 
plore other stellar contributions in §6.1). For each model, 
we obtain our full constrained luminosity function by it- 
cratively solving for the value of L*(z) that is consistent 
with our adopted ionizing rate s®(z). Of course we change 
(f)(Z) accordingly to match the SDSS normalization Eq. 2). 

The change in break luminosity for each model is illus- 
trated in Figure 7, where we plot our results at z — 3.5. 
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Model 


L 


fesc 


e LJ {z = 3,4,5) 


M*(z = 3,4,5) 


Nhdf 


NKeck 


A 


M&M-E 


0.16 


0.74, 0.40, 0.24 


-25.8, -24.6, -23.4 


U - 1O -0.07 


70 K+3V.4 
' y '°-35.8 


B 


M&M-E 


0.0 


1.44, 1.15, 0.86 


-24.2, -22.3, -20.8 


u - '°-0.29 


001+9^.4 

zoi -67.0 


C 


Proximity 


0.0 


8.25, 6.92, 2.36 


-20.6, -18.7, -18.7 


4.i8j5:S 


1790^ 



Table 1 - Our three assumed (model) cases for the ionizing background. The first column lists the model name, the second 
lists the adopted measurement of the ionizing rate, and the third column lists the escape fraction of ionizing radiation 
from stars (/ esc = implies a negligible stellar contribution to the background). The fourth column gives the weighted 
emissivity (Eq. 15) of AGN at z = 3,4,5. The fifth gives the break magnitude (Mg) for redshifts 3, 4, 5. The sixth and 
seventh column give the number of counts that would have been expected for the HDF faint AGN search (§5.2.1) and the 
Steidel search (§5.2.2) respectively. 
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Fig. 7. — The Luminosity Function for our different models at 
2 = 3.5. Notice that all three models agree at the bright-end and 
only diverge due to the location of the break. For comparison, we 
also plot the LF from MHR. 

(For comparison, we also show the MHR LF at the same 
redshift.) Note that in the limit that L m [ n <C L* <C 
L(Mb), our constraints follow the analytic relation 

log (h&)\ ~ 12.2 - 1.72 log e Q (z) - 0.81(2 - 3) 

V L ©,B / 

log ( J^=^ j ~ 2 - 80 + 2.72 log e Q (z) + 0.81(z - 3)(18) 

In practice we use these expressions as a starting point in 
our iterative solution for L* (and the implied value of <f>* ) . 

A direct presentation of each Model LF at several dis- 
crete redshifts is shown in Figure 8. For the ionizing back- 
ground of Model A, the faint-end of the LF is limited to 
virtually no evolution, while Model C allows considerable 
variation in the faint-end. The way in which these models 
evolve in terms of integrated number counts is illustrated 
in Fig. 9: for bright limiting magnitudes they show very 
similar evolution (by construction), while at faint mag- 
nitudes the evolution varies dramatically. Note that we 
have applied our emissivity constraints only for z > 3. 
For z < 2.3 we assume that the LF is fully defined by 



the 2dF results. Between these redshifts we have simply 
interpolated the values of L*, and 75. The awkward 
line shapes from z = 2.3 to 3 in Fig. 9 are due to this 
crude interpolation and should not be regarded as explicit 
predictions. 

With our Model luminosity functions in hand, we are 
now in position to discuss our results in the context of 
several faint AGN searches. Before going on to do so, we 
briefly explore the degree to which our results depend on 
our input parameter assumptions. 

5.1. Variations on Input Parameters 

The evolution of and L* for our three ionizing back- 
ground Models is illustrated in Figure 10, where we show 
our results using the parameter space introduced in con- 
junction with Figure 1 (§1). The solid, dashed, and short- 
dashed lines show how models A, B, and C evolve in this 
parameter space from z = to 6, with triangles, squares, 
and pentagons marking specific redshifts: 3, 4, and 5, re- 
spectively. All of the models overlap on the diagonal line 
representing the observed 2dF LF evolution over the range 
< z < 2.3. 

In order to estimate to what degree our input assump- 
tions might affect our results, we have recalculated our 
constraints for a series of cases in which each of our inputs 
was varied (while all other parameters were kept fixed). 
The changes we explored were: varying 7^ by ±0.2; vary- 
ing r\ v by ±40%; letting 7^ = 2.87 (rather than our fidu- 
cial 7j, = 2.58); allowing for non-zero cloud reemission of 
e R = OAi® (see §4.2); and incorporating an IGM model 
with more absorption at high-redshift (as advocated by 
Fardal et al. (1998), see Appendix). We used this suite of 
cases to estimate plausible uncertainties in our model LFs 
at each redshift, and utilize these uncertainties in the next 
subsections. The error bars on the different redshift points 
in Figure 10 represent the largest changes in </>* and L» we 
observed at each redshift for the changes described above. 
As can be seen from the figure, the primary changes are 
due to the adopted Model backgrounds, with variations on 
input parameters giving less important, although notice- 
able changes in the characteristic LF scales. 

Probably the most straightforward change we observe 
comes from from varying n u . Increasing (decreasing) r\ v 
requires fewer (more) faint AGN to reproduce the same 
ionizing background. Obtaining fewer AGN requires an 
increase in with a corresponding decrease in be- 
cause the two parameters are constrained by the SDSS 
normalization: </)» oc Ll~ Jb (recall Fig. 1). For our ±40% 
variation in rj v , the effect on <j>* is about a factor of 2 , while 
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Fig. 8. — The Luminosity Function for different redshifts in our Models A, B, and C. The lines for z = 0.35 and 2.3 are set by the 2dF data 
and, therefore, are common to all three models. 



L* varies by about 50%. Adding IGM reemission acts in 
the same way as increasing the value of r\ v . In going from 
no reemission to e R — 0A£® , increases by a factor of 
~ 2, and (f>* decreases by a factor of ~ 3. Employing a 
model with more attenuation at high redshift (Fardal et 
al. 1998, Appendix) requires more AGN to match a given 
ionizing rate. The effect is most prominent in Model A for 
z > 4. It increases 0* by nearly a factor of 6. 

We have no observational constraints on 7/ at high red- 
shift, so variations in this parameter are certainly impor- 
tant to explore. A simple approximation of the emissivity 
integral, along with the SDSS normalization, gives us the 
variation to our fiducial (7/ = 1.58) break magnitude: 

M* B (V) ~ A£(1.58) » 1.6 - 4.3 log ( ^fj^ ) • (19) 

By iterating over the relevant equations, we obtain more 
precise estimates, which we list in Table 2 for jf = 1.38, 1.78 
in each of our models. A flatter (steeper) faint-end slope 
requires more (less) AGN at "medium" luminosities to 
match a given background. This alteration is one of the 
dominant variations on Model A, since the break here is 
very near to the limits from the SDSS constraint (Eq. 2), 
and thereby the corrections to Eq. 19 are somewhat larger. 
In fact, at z ~ 3, a steeper 7/ would make M B < —26. A 
break at such a large absolute magnitude would presum- 
ably have been detected by WH094 and SSG. Indeed, if 
current indications are correct and M B > —26, then our 
fiducial Model A is close to being a minimum allowable 
LF. This possibility will be examined further in §6. 

One last parameter that we give special attention to is 
7b . The SDSS value of 2.58 is quite a bit smaller than 
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redshift 

Fig. 9. — The evolution of the integrated number counts for the 
separate models, with lines defined the same as in Fig. 7. 



in previous determinations. Similar to the case of 7/ for 
Model A, taking the steeper value of 2.87 from SSG has a 
relatively large effect on Models B and C (recall Fig. 1). 
In Model C the effect is nearly an order of magnitude in 
L*. We plot this in Fig. 10 with a separate (dotted) error 
bar, although it is only visible above the (solid) error bar 
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Model 


7/ 


M* (■/ — 3 4 *i1 


N HDF 


AT ^ , 

1 v Keck 


A 


1.38 


-25 2 -24 1 -22 Q 


0.15 


66.0 


1.78 


-27.0, -25.7, -24.3 


0.17 


90.6 


B 


1.38 


-23.7, -21.9, -20.4 


0.79 


219 


1.78 


-25.2, -23.1, -21.4 


0.67 


231 


C 


1.38 


-20.3, -18.4, -18.4 


4.64 


1990 


1.78 


-21.2, -19.2, -19.2 


3.55 


1520 



Table 2 - We recalculate each model with a higher and lower value of the faint-end slope. In comparison to our fiducial 
values in Table 1, the changes to the break magnitude are well-approximated by Eq. 19. Notice, however, that the faint 
counts (see §5.2.1 and §5.2.2) are not highly-dependent on the choice of 7/. 
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log(L./L s ) 

Fig. 10. — Solutions to the emissivity constrained models for < 
z < 6. The straight diagonal line is from the 2dF LF (0 < z < 2.3), 
which we force all of our models to agree with. The solid error bars 
come from varying parameters in the theory (see text). The dotted 
line error bars are specifically from taking the SSG value for 75. 

in cases where it exceeds all other uncertainties. Note that 
since the median redshifts of SSG and SDSS are, respec- 
tively, ~ 3.3 and ~ 4, it may be that is evolving over 
these epochs. 

5.2. Faint Surveys 

There are several observations (either planned or com- 
pleted) that are relevant to faint magnitude AGN. The 
counts obtained for the deeper, smaller fields that have 
been done so far are not as statistically significant as those 
based on large, shallow surveys we used to normalize our 
Model LF's at the bright-end. However these fainter sur- 
veys do provide a useful consistency test for our Models, 
and forthcoming surveys for faint AGN will offer further 
tests of our input assumptions about the ionizing back- 
ground. In this section we briefly investigate several sur- 
veys in light of our Model expectations and also predict 
what AGN counts might be seen by future observational 
programs. 

5.2.1. Hubble Deep Field 



The Hubble Deep Field is approximately 2.3 by 2.3 arcmin 2 . 
Looking only at unresolved sources, Conti et al. (1999) 
claim that there are no z > 3.5 QSO candidates in the 
HDF down to estimated completeness limits of B450 ~ 
27,0, V 6 o 6 - 27.0, I 8 i4 ^ 26.0. 10 The Conti et al. results 
are consistent with a similar HDF search by Elson, San- 
tiago, & Gilmore (1996). Conservatively, we will assume 
that this implies an upper bound of 3 AGN per HDF field. 
For Poisson statistics, this corresponds to 5% probability 
of seeing zero AGN. At slightly fainter magnitudes, Jarvis 
& MacAlpine (1998) take resolved sources in the HDF and 
select those with nuclei having the expected AGN colors. 
In the magnitude range V606 ~ 27.0 — 28.5, for z > 3.5, 
they find 12 candidates, but they make it clear this is an 
upper limit. We illustrate the observational bounds along 
with our three model expectations in Figure 11. Interest- 
ingly, there is one high redshift (z w 3.5) QSO seen by 
Chandra in the HDF (CXOHDFN J123639.5+621230.2, 
see Hornschemeier et al. 2001; Brandt et al. 2001). This 
object is not selected by either Jarvis & MacAlpine (1998) 
or Conti et al. (1999), even though it has R = 24.3. 

Using the Conti et al. limits (V606 < 27.0, z > 3.5), 
we have calculated what each of our Models would have 
expected in the same field. These values are listed under 
Nhdf in Table 1, along with an uncertainty calculated 
by the same variation of input parameters discussed in 
the previous section. Model C predicts too many QSOs. 
Even allowing 7& to be flatter than the SDSS value cannot 
reduce the expected value below ~ 3 per field. We discuss 
the implications of this in §6. 

5.2.2. Keck QSOs 

Unlike most other surveys, which have targeted QSOs by 
their strong UV emission, Steidel and collaborators have 
only taken spectra of objects with strong continuum breaks 
at the rest-frame Lyman limit. Their technique was de- 
signed to find high-z galaxies, which have an intrinsically 
softer spectrum at the Lyman continuum. Type I QSOs 
have much harder continua, so they would be generically 
missed by the Lyman-break selection technique. How- 
ever an estimated ~ 60% of the lines of sight to z ~ 3 
QSOs should have a foreground Lyman limit system at 
high enough redshift to create an "artificial" break. These 
are the AGN that get targeted for spectroscopic identi- 

10 Besides the attenuation from continuum absorption, the observed 
spectrum of a distant QSO will be affected by line blanketing in the 
IGM, for which we employ the approximations from Madau (1995). 
The resulting QSO colors match qualitatively the color-color plots 
from Conti et al. (1999), Haiman, Madau, & Loeb (1999), and Jarvis 
& MacAlpine (1998). 
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Fig. 11. — Number counts in the Hubble Deep Field for z > 3.5. 
Shown are the predictions from our three models, as well as the 
upper limits from two searches for QSO candidates in the HDF. 




Fig. 12. — Predicted number counts in the 1 — 3.5^tm band for 
z > 5 for our three models, as well as for MHR. The vertical lines 
are the proposed flux limits for future space telescopes. 



fication. Steidel et al. (2002) identified 13 type I QSOs 
between redshifts of 2.7 and 3.3 down to 1Z — 25.5. Based 
on the size of their field and a rough estimate of their in- 
completeness, the implied number density of type I AGN 
is approximately 190 per square degree. 1 

We have compiled the expected number per square de- 
gree in each of our Models in Table 1 (with uncertainties 
determined as before). Again, Model C seems to over- 
predict (dramatically) the number of high-z QSOs, while 
Model B, is consistent, at least within reasonable uncer- 
tainty. Interestingly, the LF motivated by the Steidel 
group's own estimate for the stellar escape fraction (Model 
A) seems to under predict the expected QSO abundance 
by a factor of ~ 2. We discuss this in the next section. 

Also detected in the same survey were 16 narrow lined, 
type II QSOs. These locally absorbed objects are quite 
important to the X-ray and, perhaps, IR backgrounds, but 
they are likely very weak emitters at ionizing frequencies, 
and therefore do not directly relate to the constraints we 
have set up. We will also address this population in §6. 

5.2.3. Future Surveys 

There are a number of future and ongoing surveys that 
plan to fill in the faint high redshift portion of the AGN 
LF. The Big Faint Quasar Survey (Hall 2001) will presum- 
ably probe the LF for -26.5 < M B < -23.5, for z > 3.5. 
This is below the luminosity range of SDSS. Hall predicts 
that this survey (for B < 23.5 and 7deg 2 ) will find - 120 
quasars for 3.5 < z < 4, and ~ 40 for 4 < z < 5. While 
Model A with a substantial stellar contribution would pre- 
dict fewer detections than this, Models B and C with AGN- 
dominated backgrounds expect 2-3 times as many QSOs. 

The Space Infra-Red Telescope Facility (SIRTF) is pro- 
posed to have l^Jy sensitivity, and the Next Generation 

11 Steidel et al. (2002) state that more careful analysis will be done 
in a later paper. 



Space Telescope (NGST) is proposed to have InJy sensi- 
tivity between 1 and 3.5/im. Although selecting high red- 
shift QSOs in the Infra-Red may not be straightforward 
(see Warren & Hewett 2002), it is worthwhile to work out 
our theoretical predictions at these detection thresholds. 
These are shown in Figure 12. 

Since all of our Models drop off steadily for z > 5, the 
number counts illustrated in Figure 12 are dominated by 
the z ~ 5 QSOs. At this redshift for our chosen SED 
and cosmology, the proposed flux limits of SIRTF (NGST) 
would correspond to AGN with L ~ 10 n L o (10 8 L©). We 
can compare our results to those of Haiman & Loeb (1998), 
who used a model based on CDM structure formation to 
predict that SIRTF would observe on the order of a few 
z > 5 AGN, while NGST would see approximately 10 3,5 . 
A close examination of Haiman & Loeb's differential LF 
(their Fig. 5) at L ~ 10 11 Lq^b shows that it is an order of 
magnitude higher than the simplest extrapolation of the 
SDSS results. This explains why their SIRTF prediction 
is more than an order of magnitude higher than any of our 
predictions. So too, their NGST estimates are a factor of 
ten to a hundred above any of our estimates. However, this 
appears to be due to the fact that their LF is practically 
constant for L < 1O 1O L and z > 5. 12 

6. IMPLICATIONS 

In addition to direct constraints on the nature of the 
faint- AGN LF, our results have some interesting implica- 

12 Note that Haiman &; Loeb have since re-tooled their model 
(Haiman, Madau & Loeb 1999; Haiman & Loeb 1999), further sup- 
pressing the formation of small mass black holes, and the updated 
NGST predictions appear to be more in line with our Models B and 
C. Still, however, the predicted shape of the bright end of the LF 
from a simple scaling of the CDM mass function would appear to 
be in conflict with the SDSS observations, and therefore, it may be 
necessary to employ some sort of feedback mechanism to curtail the 
formation of intermediate AGN. 
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tions for the nature of ionizing sources as well as some 
typical characteristics of AGN themselves. 

6.1. Ionizing Sources 

Our joint analysis of the ionizing rate and faint AGN 
counts provides a potentially useful avenue for exploring 
the contribution of stars and AGN to the ionizing back- 
ground. We showed in §5 that faint surveys do not find 
the number of QSOs that would be expected if all of the 
proximity-effect-derived background is coming only from 
AGN (Model C). Conversely, a model that reproduces the 
(lower) ionizing background intensity favored by the flux 
decrement technique under-predicts the faint counts un- 
less the stellar contribution is limited (Model B). In this 
subsection we attempt to place more quantitative limits 
on the stellar contribution implied by the two competing 
background intensity measurements. 

For each ionizing background measurement, we have cal- 
culated the expected number of Keck and HDF counts for 
several values of the escape fraction (§4.1). The results 
are shown in Figure 13, with the different symbol types 
reflecting different background assumptions: the squares 
represent a standard proximity effect background (Fardal 
et al. 1998) and the circles correspond to an assumed 
flux decrement background. In order to be conservative, 
we also include results for a background that matches the 
lower limit on the proximity effect rate taken from Scott 
et al. (2002) (triangles), with T_i2 ~ 1.0 over the range 
1-7 % z < 4.1. Solid symbols include no IGM reemission 
and open symbols include 40% cloud reemission. The hor- 
izontal dashed lines indicate the upper limit from the HDF 
search (upper panel) and observed AGN per square degree 
seen in the Keck fields (lower panel). 

By examining Figure 13 we see that even with (the un- 
physically large value) / csc = 1, the standard proximity 
effect background yields far too many Keck QSOs. A con- 
tribution from IGM reemission (§4.2) does little to lessen 
the discrepancy. The lower bound on the proximity ef- 
fect from Scott et al. provides more reasonable numbers, 
but does require a significant stellar component, / csc a20%, 
even if we allow for a large amount of reemission. 13 The 
flux-decrement-derived background, on the other hand, is 
not compatible with much stellar contribution: / csc s5% is 
required in order to reproduce the Keck counts. 14 

An additional way to examine the question of stellar es- 
cape fraction involves using the He Lya forest to limit the 
relative amount of ionizing radiation from stars and AGN. 
Stars, unlike AGN, do not have much emission that ex- 
tends to the helium Lyman continuum. This means that 
a stellar-dominated background should have a smaller rel- 
ative fraction of fully ionized helium. At z ~ 2.5, the 
implied optical depth to Hen absorption is low enough 
that the background at these redshifts should be QSO- 
dominated (Davidson et al. 1996). But at slightly higher 
redshifts, Heap et al. (2000) appear to detect a Gunn- 
Peterson trough blueward of helium Lya in the spectra 
of Q0302-003. They conclude that the ratio of hydro- 

13 A similar discussion on acceptable escape fractions in light of 
proximity effect measurements can be found in Bianchi et al. (2001), 
but there they assume PLE for their AGN LF, similar to that of 
MHR. 

14 We stress that these escape fraction limits cannot be more precise, 
as the uncertainties in the Keck survey are not yet known. 
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Fig. 13. — The expected number counts in the HDF (top) and 
Keck (bottom) surveys as a function of the escape fraction. The 
upper bound on the HDF and the observed number from Steidel et 
al. (2002) are plotted as dashed horizontal lines. The squares are for 
an ionizing background from the proximity effect (as parameterized 
by Fardal et al. 1998); the triangles are from the lower bound on the 
proximity effect from Scott et al. (2000); and the circles are from the 
flux decrement measurement of M&M-E. The open/solid points are 
with/without the contribution from cloud reemission (i R = OAs®). 
Note that we do not plot points at high escape fraction for the 
flux decrement (and for the lower bound on the proximity effect, as 
well), since the resulting LF breaks would contradict observations 
(see §5.1). 



gen to helium ionization rates rises abruptly at z w 3 to 
rm/rHoii — 800, suggesting that a soft stellar contribu- 
tion is beginning to dominate the hard ionizing spectrum 
of AGN at this epoch. This kind of rapid hardening of 
the background at z ~ 3 is supported by the analysis of 
Songaila (1998) of the Siiv/Civ ratio in the IGM. Al- 
though we do not model the propagation of helium ioniz- 
ing photons in a very sophisticated way (see Appendix), we 
show in Figure 14 that Model A reproduces the reported 
evolution of rHi/rfjeii in a qualitative way. A more careful 
analysis of this phenomena will likely have to include the 
concomitant reionization of Hen (see Theuns et al. 2001). 
Recently, Sokasian, Abel, &. Hernquist (2002) have used 
observed opacities of HI and Hell to conclude that stars 
and AGN must contribute roughly equally to the ionizing 
background at z ~ 3. Theuns et al. (2002) and Bernardi 
et al. (2002) argue that the evolution of the flux decre- 
ment distribution studied via a sample of SDSS quasars 
points to Hell reionization at z ~ 3 — 4. 15 

15 In writing this paper, we explored another idea for separating the 
relative galaxy/AGN contribution to the background. This was to 
use ~ GeV gamma ray attenuation to measure the UVB as is done 
for the IR background using ~ TeV sources (see Primack et al. 2001, 
Bullock et al. 2002). Unfortunately, at these redshifts, UV photons 
are completely overwhelmed by foreground optical photons, because 
the peak energy of the interaction scales as (l+z)~ 2 . One can phrase 
this result in the positive: gamma ray attenuation is sensitive only 
to the integrated stellar light over the history of the universe. AGN 
consistent with current estimates of the ionizing background will 
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Fig. 14. — The softness parameter, Fm/Tilen, for Model B (QSOs 
only) and Model A (stars and QSOs). The rapid rise at z ~ 3 for 
Model A is qualitatively similar to that reported by Heap et al. 
(2000). 

6.2. AGN Characteristics 

In addition to quantifying our expectations for AGN 
counts in the presence of different ionizing backgrounds, 
our results have implications for the emission efficiency 
and lifetimes of AGN. 

As discussed in §1, AGN activity is likely driven by ac- 
cretion onto super-massive black holes. Under this as- 
sumption, the present-day density of black holes is re- 
lated to the integrated AGN emissivity over the history 
of the universe, modulo the AGN efficiency (see Soltan 
1982; Chokshi & Turner 1992; Haehnelt, Natarajan & Rees 
1998). Because our models provide estimates of the max- 
imum AGN emissivity out to high redshifts, we can con- 
volve our results with an estimate of the density of black 
holes today in order to determine an implied efficiency. 

Specifically the relic black hole density p. is related to 
the integrated QSO energy density by pq = ep., where 
e = Lboi/Mc 2 is the QSO efficiency of converting mass 
into energy. The total energy output of QSOs can be 
obtained by integrating the emissivity over the history 
of the universe. If we define £b(z) — vb£b{vb, z) and 
vb = c/4400A, then we can write 

PQ = — FT — T f -j-dz ■ eb{z). (20) 

Because we are interested in the total energy output, we 
have applied a bolometric correction: L^i/ 1 vbL v {vb) = 
11.8 (Elvis et al. 1994). We plot the B band emissivity 
for our three models in Figure 15. For z < 2.3 we have 
used the 2dF LF fit from Boyle et al. (2000) to determine 
the emissivity, and between z = 2.3 and 3, we linearly in- 
terpolate from the 2dF fit to our high redshift constraints. 

provide a negligible contribution to the extra-galactic background 
light relevant for gamma-ray attenuation. 



When we integrate from z = to 6, we find that Models 
A, B, and C give pq = 1.38, 1.57, and 2.81 respectively, in 
units of 10 4 M Q Mpc" 3 . The low-z LF from the 2dF alone 
gives pq(z < 2.3) = 1.12 in the same units. Thus a large 
fraction of the energy output from AGN seems to have oc- 
curred at late times even for our most extreme Model (C), 
and more than half for Models A and B . 

The present day density in black holes can be estimated 
by assuming a typical black hole mass to spheroid mass 
ratio, M,/M sp , and combining it with a local determina- 
tion for the mass fraction in spheroids f2 sp . Adopting the 
fiducial value of M./M sp = 0.13% (Merritt & Ferrarese 
2001; McLure & Dunlop 2002; van der Marel 1999) and 
tt sp = 0.002/i" 1 (Fukugita et al. 1998), we obtain 

(21) 

Salucci et al. (1999) obtained a similar value (8.2h 2 x 
10 5 M Q Mpc~ 3 ) by convolving a distribution of M,/M sp 
values with an estimate of the spheroid mass function 
(see also Yu & Tremaine 2002). Adopting p. — 5.0 x 
10 5 M Q Mpc- 3 , we obtain e = p Q /p. = 0.028,0.032, and 
0.056 for models A, B, and C respectively. 

Note that although these are in principle maximum effi- 
ciencies, based on maximum allowable AGN emissivities 
(because we have ignored reemission), they are signifi- 
cantly smaller than the often-adopted value of ~ 0.1. In 
addition, our derived efficiencies are significantly smaller 
than those obtained using the hard X-ray emissivity (Fabian 
& Iwasawa 1999; Salucci et al. 1999; Elvis et al. 2002). 
One possible explanation for this discrepancy is the exis- 
tence of a large population of obscured AGN. If AGN ef- 
ficiencies are typically zl0% (e.g. Elvis et al. 2002), then 
our results require that z 50 — 70% of AGN are signifi- 
cantly obscured in the UV and optical. Interestingly, this 
result is similar to that obtained by synthesis modelers, 
~ 75 — 90%, based on the shape of the X-ray background 
spectrum (Comastri et al. 1995; Gilli et al. 2000) . 

We have not considered type 2 AGN in our analysis sim- 
ply because the implied absorption would mean they are 
not strong ionizing sources. But unification models as- 
sume that the type 2 population is not separate from the 
type Is, but merely the result of an orientation effect. This 
would imply a non-trivial mapping between the mass func- 
tion of supermassive black holes and the optical luminosity 
function. Therefore, to study the accretion history of the 
universe, it would seem to be more straightforward to use 
bands less affected by obscuration (see Marconi & Salvati 
2001; Barger et al. 2001). Unfortunately, the X-ray LF at 
high redshift is limited statistically (Miyaji et al. 2000), 
and the AGN activity in many IR sources is still a matter 
of debate (see Lawrence 2001; Fadda et al. 2002). At high 
redshift, the best limits on the shape and evolution of the 
AGN LF still come from the optical/UV. 

Another constraint related to the local number density 
of black holes concerns the fraction of supermassive black 
holes that are active at any time, f oa (z). This can be 
related to a typical AGN lifetime, t agn , which we define, 
via f n(z) = t agn {z)/tji{z), to be the time that a typical 
black hole is active over the course of a Hubble time. For 
simplicity, let us assume that every AGN shines at a fixed 
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Fig. 15. — The cmissivity in the B band. At low redshift we simply 
take the fit from 2dF. For z > 3, we plot our three models with a 
simple linear interpolation to the 2dF limit of z = 2.3. We have 
adopted units such that one Lq b corresponds to 2.11 X 10 33 erg/s 
(see SSG). 



fraction of the Eddington luminosity: L — L\ = AL^dd, 
where L E dd ~ 6 x lO 3 (Af./M )L Qi B- With this assump- 
tion, we write / on (averaged over luminosity) as 



/on (z) 



$(> L\, z) 



> 



*(> L x ,z) 



(22) 



N(> M.,z) ~ N{> M.,z = 0) 
The inequality comes from the assumption that a given 
black hole's mass only grows with time and that the black 
hole mass function is a decreasing function of mass. One 
possible caveat to this inequality is that black hole merging 
could (but not necessarily) reduce the number of low-mass 
black holes with time. In this case, the above inequality 
would break down. However, in the case of halos of mass 
~ 10 10 M© (corresponding to black holes of ~ 10 6 M Q ), 
the number density does not decrease substantially from 
z ~ 5 to the present. This is because most merging occurs 
in high-mass-ratio events. 

So assuming Expression 22 is valid, we can take for the 
denominator the results of Salucci et al. (1999), who es- 
timate that N(> M.,z = 0) ~ lCT 2 Mpc~ 3 for M. = 
10 6 M Q . Given A, this number allows us to determine 
the implied lower limits on f on {z) for each of our Mod- 
els. We plot these limits in Figure 16 for A = 1.0 (L\ = 
6 x 1O 9 L 03 ) and A = 0.05 (L x = 3 x 1O 8 L , B ). Fortu- 
nately, none of our models result in the unphysical / on > 1. 
But for the corresponding lifetime limits, Model C requires 
very long-lived AGN at high-z: t agn al0 8 yr, whereas Mod- 
els A and B suggest lifetimes that are on the high-side of 
what is typically assumed: t agn ^10 7 yr. If estimates of the 
AGN lifetime from other lines of reasoning can be obtained 
(e.g. Martini & Weinberg 2001; Haiman & Hui 2001), then 
this sort of analysis may be capable of restricting theories 
on the history of black hole accretion (see Ciotti et al. 
2001). 
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Fig. 16. — Top- lower bounds on the fraction of 10 6 Mq black 
holes that are active at any redshift; bottom- implied lower bounds 
on typical AGN lifetimes (see text for a description). We plot two 
assumptions for the ratio of AGN luminosities to the Eddington 
luminosity: A = hj L^dd equal to 0.05 (upper set of curves) and 1.0 
(lower). Our three models are represented, along with the model of 
MHR for comparison. The lines are defined as in Figure 7. 



7. SUMMARY AND CONCLUSIONS 

In this paper, we have explored the implied evolution 
of the faint AGN luminosity function using three different 
models for the ionizing background from 3 &z s6. Our LFs 
are derived by matching the implied faint-end emissivity. 
Unfortunately, the value of the ionizing background rate 
is not universally agreed upon. Measurements obtained 
using the proximity effect generally give larger values than 
those based on the flux decrement distribution. 

Although the current data on faint AGN come from 
small-field searches with only modest statistical accuracy, 
we have used them to obtain rough evaluations of our Mod- 
els. If AGN were producing a background at the level mea- 
sured by the proximity effect technique, then significantly 
more would have been seen in the HDF and LBG fields. 
For typical proximity effect values of the background, even 
an escape fraction of unity would require more than three 
times the number of AGN observed in the Keck fields. We 
were able to obtain modest agreement only by taking the 
lowest bounds on the proximity effect measurements and 
by including a high galactic escape fraction: / csc ft20%. 
Conversely, if the flux decrement rate is adopted, there is 
little room for a significant stellar component, and / csc s5% 
is required to match the faint AGN counts. Future AGN 
searches and developments in our ability to measure the 
background intensity will be useful for further constrain- 
ing the stellar contribution to the ionizing background at 
high-z. 

We used our derived luminosity functions and local de- 
terminations of the black hole relic density to determine 
the typical AGN efficiency of converting mass to light. Al- 
though many previous estimates based on X-ray counts 
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have obtained e zO.l (e.g. Elvis et al. 2002), we find sig- 
nificantly lower values e ~ 0.028 — 0.056, suggesting that 
more than half of all AGN are obscured in the UV/optical. 
Alternatively, lower than expected values of e may be de- 
rived if much of a black hole's mass is set by a massive 
"seed" before the onset of its active phase. A joint analy- 
sis of UV-background measurements, optical counts, and 
X-ray observations would likely resolve this issue. 

A similar analysis allows us to limit the fraction of black 
holes that could have been active at any redshift. For our 
models that are most consistent with faint AGN counts, 
we get / n £0.5% at z ~ 2.5. We can interpret this in terms 
of a lower limit on the average AGN lifetime: t agn £l0 7 yrs. 
For comparison, an AGN luminosity function that matches 
the proximity effect background (but over-produces the 
faint counts) requires / on -5% and t agn £l0 8 yrs. Note that 
these numbers are firm lower limits because we assume 
that there are no obscured AGN. If a fraction of AGN 
were obscured in the optical, then these limits would be 
increased by the inverse of the same fraction. 

Finally, we conclude with some remarks about CDM- 
based theories of AGN formation and how they might 
be constrained. Usually cosmological models of this kind 
rely on simple mappings between dark halo masses and 
AGN luminosities, but a fundamental unknown concerns 
the amount of feedback that takes place (e.g. Haehnelt & 
Rees 1993, Haiman & Loeb 1998). We point out that one 
potentially worrisome difficulty for models with very little 
feedback arises from direct observations of the AGN LF 
at bright luminosities. Specifically, the bright-end slope 
measured by the SDSS at high redshift is even flatter than 
that measured by the 2dF at low redshift. This is the op- 
posite of what is expected if the AGN LF maps simply to 
the halo mass accretion rate function (see, e.g., Haiman & 
Loeb 1998). Therefore the SDSS+2dF results alone seem 
to suggest that some feedback is needed for bright to in- 
termediate luminosity AGN. 

Another question is whether any (additional) feedback 
might be needed in low-mass/ low-luminosity objects in 
order to explain the observed ionizing background rate at 
high-z. As discussed above, our Models A and B span 
what we believe to be the range allowed by the ioniz- 
ing rate measurements. Model B corresponds to the case 
where AGN dominate the background, and, interestingly, 
its evolution is very reminiscent to that of the halo mass 
function in CDM: faint AGN are numerous at early times, 
but their numbers fall-off slowly at late times (see Figure 
8). It would seem that if AGN do dominate the ioniz- 
ing background, then CDM-models would require no addi- 
tional feedback on low mass scales (relative to what might 
already be required at the high-mass end in order to match 
the SDSS results.) 

If instead stars contribute with an escape fraction con- 
sistent with the Steidel et al. (2002) report (Model A), 
then many fewer AGN are permitted, and CDM-based 
models would need to be adjusted further in order to ad- 
ditionally suppress the formation of low-luminosity AGN. 
Both Haiman, Madau & Loeb (1999) and Kauffmann & 
Haehnelt (2000) presented models with this kind of luminosity- 
dependent feedback, but it is unclear if, as presented, these 
models would be able to reproduce both the ionizing rate 
(with an allowance for stars) and the faint AGN counts in 



the HDF and LBG fields. It would be interesting to see 
such a comparison. The model of Kauffmann & Haehnelt 
(2000) is ideally suited for this test because it predicts the 
joint population of galaxies and AGN in a self-consistent 
manner. 16 

It is becoming increasingly clear that galaxy formation, 
AGN activity, and the resulting ionizing background are 
intimately connected and take part in an important feed- 
back loop (Gebhardt et al. 2000a; Ferrarese & Merritt 
2000; Kauffmann & Haehnelt 2000; Bullock, Kravtsov, & 
Weinberg 2000; Benson et al. 2002; Somerville 2002). For 
this reason the need is high for self-consistent models that 
treat all of these processes together, as is the need for ad- 
ditional observational constraints on models of this kind. 
The limits presented here provide one small step in this di- 
rection, but if some agreement can be reached on measure- 
ments of the ionizing background emissivity, the resulting 
constraints would prove remarkably important for piecing 
together the story of cosmological structure formation. 
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are less abundant and therefore are less important to the 
effective attenuation. In any case, a more careful study of 
the helium ionization state should take into account the 
self-shielding and reemission in the IGM. 



APPENDIX 

The following comes primarily from MHR, but see also 
Fardal et al. (1998). The attenuation due to photoelectric 
absorption can be modelled by a distribution of discrete 
absorbers (or "clouds"): 

r efM ^) = l i*l dAtajj^d-.-^) (1) 

where v' = v (1 + z')/(l + z), and, for a given absorber, 
A^hi is the column density in neutral hydrogen. This dis- 
tribution is usually parameterized as: 

d 2 N 1 fN w x ~ 



- ) (1 + ^) 7 (2) 



If we assume that the distribution is smooth over red- 
shift and column density, then we get an analytical ex- 
pression for the effective optical depth. We assume N = 
1.6 x 10 15 cm- 2 , (3 = 1.5, 7 = 2, which MHR claim ap- 
proximate more precise formulations. [In §5.1 we compare 
what our results would be using the corresponding model 
from Fardal et al. (1998): N a = 3.27 x 10 14 cnr 2 , = 1.5, 
7 = 2.58.] As for the attenuation in a particular cloud, we 
will consider only Hi and Hen absorption: 

T(v') = N m <T m {v') + ^HcIIffHcIl(^) (3) 

The cross sections are zero below threshold (fjj = 13.6eV; 
t'Heii = 54.4eV), and above threshold are proportional to 
z^~ 3 , with cth(^h) = 6.35 x 10~ 18 cm 2 and <thcIi(^HcIi) = 
1.59 x 10~ 18 cm 2 . These cross sections are separated far 
enough in frequency that we can consider the ionization 
state of hydrogen and helium separately. We can make 
the approximation that the relevant clouds are optically 
thin (see Haardt & Madau 1996; Fardal et al. 1998), so 
that: 
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This relation will not hold true in high column density 
clouds because of self-shielding effects, but these clouds 



